On the Retrieval of Sea Ice Thickness and Snow Depth using Concurrent Laser Altimetry and L-Band Remote Sensing Data

The accurate knowledge of sea ice parameters, including sea ice thickness and snow depth over the sea ice cover, are key to both climate studies and data assimilation in operational forecasts. Large-scale active and passive remote sensing is the basis for the estimation of these parameters. In traditional altimetry or the retrieval of snow depth with passive microwave sensing, although the sea ice thickness and the snow depth are closely related, the retrieval of one parameter is usually carried out under assumptions over the other. For example, climatological snow depth data or as derived from reanalyses contain large 5 or unconstrained uncertainty, which result in large uncertainty in the derived sea ice thickness and volume. In this study, we explore the potential of combined retrieval of both sea ice thickness and snow depth using the concurrent active altimetry and passive microwave remote sensing of the sea ice cover. Specifically, laser altimetry and L-band passive remote sensing data are combined using two forward models: the L-band radiation model and the isostatic relationship based on buoyancy model. Since the laser altimetry usually features much higher spatial resolution than L-band data from Soil Moisture Ocean 10 Salinity (SMOS) satellite, there is potentially covariability between the observed snow freeboard by altimetry and the retrieval target of snow depth on the spatial scale of altimetry samples. Statistically significant correlation is discovered based on high-resolution observations from Operation IceBridge (OIB), and with a nonlinear fitting the covariability is incorporated in the retrieval algorithm. By using fitting parameters derived from large-scale surveys, the retrievability is greatly improved, as compared with the retrieval that assumes flat snow cover (i.e., no covariability). Verifications with OIB data show good 15 match between the observed and the retrieved parameters, including both sea ice thickness and snow depth. With detailed analysis, we show that the error of the retrieval mainly arises from the difference between the modeled and the observed (SMOS) L-band brightness temperature (TB). The narrow swath and the limited coverage of the sea ice cover by altimetry, as well the uncertainty associated with the radiation model are potential sources of error. The proposed retrieval algorithm (or methodology) can be applied to the basin-scale retrieval of sea ice thickness and snow depth, using concurrent passive remote 20 sensing and active laser altimetry based on satellites such as ICESat and ICESat-2. 1 The Cryosphere Discuss., https://doi.org/10.5194/tc-2017-133 Manuscript under review for journal The Cryosphere Discussion started: 9 August 2017 c © Author(s) 2017. CC BY 4.0 License.